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Finite temperature 2-color QCD for real and imaginary chemical potential 
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In this work we study the finite temperature SU(2) gauge theory with staggered fermions for real and imaginary 
chemical potential. We test the method of analytical continuation of Monte Carlo results obtained for imaginary 
chemical potential, by comparison with those obtained directly for real chemical potential. 
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1. INTRODUCTION 

Understanding the QCD phase diagram in the 
temperature-chemical potential (T, fi) plane is 
important for the implications in cosmology, as- 
trophysics and heavy-ion collisions, where quark- 
gluon plasma could be formed. The main purpose 
of the most recent studies is to localize the crit- 
ical lines in the phase diagram and to determine 
the kind of transition across them. The most 
powerful tool for such investigation is the for- 
mulation of the theory on a space-time lattice. 
Unfortunately, at nonzero chemical potential the 
determinant of the fermion matrix becomes com- 
plex, thus requesting nonstandard approaches to 
Monte Carlo simulations. 

One way to circumvent this problem is based 
on the idea to perform numerical simulations at 
imaginary chemical potential ^ = and to 
analytically continue 2 3 4 5 the results to real ii. 
The advantage is that for imaginary fi the fermion 
determinant is again real, so simulations are fea- 
sible in QCD. Moreover, this method allows to 
check the results obtained by other approaches. 
On the other side the method is applicable as long 
as one remains inside a region of analyticity of the 
partition function in the {T,fij) plane. 



2. PHASE DIAGRAM 

Roberge and Weiss (RW) have shown P that 
the partition function of the SU(Nc) theory is 
periodic in the parameter 6 = iJ-i/T, owing to 
Z(Nc), and have observed that there are first- 



order vertical lines in the (T, 9) plane, located at 
9 = 27r(/c -I- 1/2) /Nc, with integer k. According 
to RW, these transitions are placed above a cer- 
tain temperature Te and since the first of them 
appears for 9 = tt/Nc, they conclude that there 
is not transition at /i = 0. The last conclusion is 
contradicted in the literature [7| where there are 
different arguments according to which a chiral 
phase transition for /i = does exist. Its na- 
ture depends on the number of quark, on fermion 
masses, and so on. 

In view of this, the worst case which can oc- 
cur, for the purposes of analytical continuation 
(AC), is that the RW critical fine at 9 = n/Nc 
could be prolongated downwards and bend up 
to reaching the 9 — axis in correspondence of 
a critical temperature Tc- The RW periodicity 
and the fj, — > —fj, symmetry, resulting from CP 
invariance, imply that knowing the expectation 
value of any observable for T and 9 inside the 
strip {0 <9 < tt/Nc, < T < oo} in the {T,e) 
plane is enough to fix it on the remaining part of 
the plane. It is more suitable for lattice simula- 
tions to introduce the inverse coupling P and the 
chemical potential in lattice units, (1 = a/i. The 
RW transition lines in the {(3,fii) plane appear 
at /i/ = 2'n{k -1-1/2)/ (NcNr) and we can consider 
only the strip {0 < /i/ < n/{NcNr); < /3 < oo} 
(see Fig. Pi . 

3. NUMERICAL RESULT 

In this work we consider SU(2) gauge theory 
with staggered fermions {rif =8) because in this 
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Figure 1. Schematization of the phase diagram in 
the (/3, jli) plane. The curved dotted hne repre- 
sents the critical line in the (/3, flu) plane. 
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case numerical simulations are feasible for both 
real and imaginary chemical potential, thus al- 
lowing to test the method of AC by direct com- 
parison. We performed all simulations on a 8'^ x 4 
lattice, setting the quark bare mass to am — 0.07. 
For different values of (3 we performed simulations 
for both real and imaginary chemical potential, 
varying /</ and jjn between and t:/{2Nt) and 
determined expectation values of the Polyakov 
loop L and the chiral condensate (ipip). Then, we 
interpolated data obtained for imaginary fi with a 
truncated Taylor series of the form a + b jij + cjlj. 
After "rotating" this polynomial to real chemi- 
cal potential, thus leading to a — & fij^ + cpi^, we 
compared it with the determinations obtained di- 
rectly for real ^ = ^r. 

We selected three values of /3 (1.90, 1.45, 0.90), 
each representative of one of the regions /3 > /3_b, 
I3c< 13 < Pe and (3 < j3c. 

For /? = 1.90 we interpolated the data obtained 



for < A/ < 7r/(2iV^) and found (see Fig. ^ 
that the rotated polynomial perfectly interpolates 
data obtained directly for real /i — in the same 
range. From this outcome we can conclude that 
the AC works in the region j3 > (3e for all the 
/i/'s before the first RW critical line. 

For P = 1.45 we are in the region f3c < P < Pe 
for which we expect that by varying /i/ the chiral 
critical line is crossed, while no transitions should 
be met at the same P by varying fiR. Therefore 
this time we interpolated data obtained for the 
real chemical potential < /ifl < tt/{2Nt-) and 
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Figure 2. Polyakov loop (above) and chiral con- 
densate (below) vs p. for P — 1.90. 



rotated this polynomial to its counterpart in /Ij. 
The comparison (see Fig. O with data obtained 
directly for imaginary chemical potential shows 
agreement for /i/ below ~ 0.18, which can be 
taken as an estimate of the chiral critical value 
in the (/?, /i/) plane at the given p. 

For P — 0.90 < Pc, we expect analyticity for 
all possible values of fii, while the chiral critical 
line in the (/3, /ir) plane is crossed by varying fiR 
at the given value of p. Therefore we interpo- 
lated data obtained for < /i/ < tt/{2Nt) and 
rotated this polynomial to its counterpart in fiR. 
The comparison (see Fig. 0)) with data obtained 
directly for real chemical potential shows agree- 
ment for jj,R below ~ 0.12, which can be taken 
as an estimate of the critical value in the (/3, (Ir) 
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Figure 3. Polyakov loop (above) and chiral con- 
densate (below) vs /i for /3 = 1.45. 
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Figure 4. Polyakov loop (above) and chiral con- 
densate (below) vs /i for /3 = 0.90. 



plane at the given (3. 

4. CONCLUSION 

We have presented a possible structure of the 
phase diagram in the (T, /i/) plane for 2-color 
QCD with Uf = 8 and have found that the 
method of AC works fine within the restrictions 
imposed by the presence of nonanalyticities. 
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